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We consider the phase separation in an asymmetrically charged lipid bilayer membrane consisting 
of neutral and negatively charged lipids that are in contact with in and out ionic solutions having 
different ionic strengths. The two asymmetric leaflets are coupled through electrostatic interactions. 
Based on a free-energy approach, the critical point and phase diagrams are calculated for different 
ionic strengths of the two solutions and coupling parameter. An increase of the coupling constant or 
asymmetry in the salt concentration between the in and out solutions yields a higher phase separation 
temperature because of electrostatic interactions. As a consequence, the phase-coexistence region 
increases for strong screening (small Debye length). Finally, possible three-phase coexistence regions 
in the phase diagram are predicted. 



I. INTRODUCTION 

Model membranes consisting of mixtures of saturated 
lipids, unsaturated lipids, and cholesterol have received 
considerable attention in recent years because they can 
be regarded as model systems of biological cell mem- 
branes. Below a certain temperature, the membrane 
undergoes a lateral phase separation between a liquid- 
ordered phase (rich in saturated lipid and cholesterol) 
and a liquid-disordered phase (rich in unsaturated lipid) . 
The resulting lipid domains are sometimes called "rafts" , 
and are believed to play an important role in various bio- 
cellular processes such as signal transduction and coop- 
erative membrane trafficking 

A large number of experimental studies have been car- 
ried out to further explore the consequences of domain 
formation in model membranes and their relation with 
biomembrane functioning. They include, among others, 
studies of domain morphology l2|, |1] , domain budding [1], 
growth dynamics of domains [5[, and formation of peri- 
odic lateral structures @. In particular, we note that 
lateral phase segregation was directly observed for lipids 
at the air-water interface and for giant lipid vesicles us- 
ing fluorescence microscopy Furthermore, in an 
attempt to understand these experimental findings, sev- 
eral theoretical models have been proposed 

More recently, several experimental works have been 
conducted in order to understand how electrostatic in- 
teractions affect the phase separation of model mem- 
branes composed of charged and neutral lipids [Tol . 
[T]| . As charged lipids are everpresent in biomem- 
branes, the role of electrostatic interactions is impor- 



tant also for biological cells. The phase-coexistence re- 
gion was reported to be fully suppressed in lipid bi- 
layers consisting of a three-component mixture of neu- 
tral saturated lipid DPPC (1,2-dipalmitoyl-sn-glycero- 
3-phosphocholine) , negatively charged unsaturated lipid 
DOPS (l,2-dioleoyl-sn-glycero-3-phospho-L-serine), and 
cholesterol, because of the electrostatic repulsion be- 
tween charged lipids Furthermore, the extent of 
the phase-coexistence region appears when salt (CaCl2) 
is added, an effect that can be understood by the screen- 
ing of electrostatic interactions in presence of added 
salt. In a related experimental study, it was shown that 
for mixtures of DPPC, negatively charged DOPG (L- 
a-l,2-dioleoyl-sn-glycero-3-phosphoglycerol) and choles- 
terol, the phase-separation temperature becomes higher 
as salt is added TT] . Such effects were also considered in 
several theoretical studies [12-; 14] . 

Another important aspect of biomembranes is the 
compositional asymmetry between their inner and outer 
leafiets. Collins and Keller investigated the phase be- 
havior in asymmetric lipid bilayers, and showed that the 
domain formation in each leafiet is coupled through an 
inter- leaflet interaction [l^. In some cases, it was re- 
ported that phase separation in one of the leaflets in- 
duces a phase separation in the second leaflet, while in 
other cases, the lack of phase separation in one leaflet 
suppresses domain formation in the second leaflet. The- 
oretical models based on regular solution theory [l6j or 
Landau theory [l3] attempted to take into account inter- 
leaflet interactions and explored their consequences on 
the bilayer phase diagram. 

In related works based on the Poisson-Boltzmann the- 
ory, a model describing the immiscibility transition in 
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asymmetric and chaumd membranes was proposed by 
May and coworkers |fl, [l^ . Their model took into ac- 
count an inter-leaflet electrostatic coupling between two 
rigid and charged planes modeling a bilayer membrane. 
The spinodal line, characterizing the phase coexistence, 
and the critical point were derived and depend on the 
ionic strength of the solution as well as on the inter- leaflet 
electrostatic coupling. One of the main conclusions was 
that a stronger inter-leaflet coupling enhances the lateral 
phase separation. However, it was equally assumed by 
the authors that the two aqueous solutions in contact 
with the two leaflets have the same salt concentration. 
We note that this assumption is an over-simplification 
for cellular biomembranes, where due to ionic channels 
and other active processes, the ionic strength is differ- 
ent on the two sides of the membrane (the inner- and 
extra-cellular regions). This difference in ionic strength, 
in turn, contributes to an additional gap between the sur- 
face potential on the two sides of the bilayer and plays an 
important role, for example, in neuro-transmission pro- 
cesses. 

In the present work, we extend the model of Refs. (Tsl . 
[l9l | in order to mimic asymmetric biomembranes in 
a more complete way. In Sec. HIl we introduce a 
model based on Poisson-Boltzmann theory to describe 
phase separation for asymmetric charged bilayers "sand- 
wiched" between two ionic solutions having different ionic 
strength. In Sec. IIIH the variation of the critical point 
and the entire phase diagram as a function of the ionic 
strength (Debye length) and inter-leaflet electrostatic 
coupling is explored. It is shown that the phase sepa- 
ration is enhanced as the salt concentration is increased, 
and is due to the enhanced electrostatic screening. When 
an electrostatic coupling between the leaflets is intro- 
duced, the lateral phase separation occurs at higher tem- 
peratures as compared with the no coupling case. We also 
show that the phase-separation temperature increases 
when the concentration difference between the two salt 
reservoirs becomes larger. Finally, discussion and com- 
parison with other works are presented in Sec. IIVI 



II. MODEL 

The system is modeled as a binary lipid-bilayer com- 
posed of a mixture of negatively charged and neutral 
lipids. Although many experiments are done in the pres- 
ence of added cholesterol as a third component, it is rea- 
sonable to stay within the simpler case of a binary lipid 
mixture. As cholesterol is non-charged, its presence will 
not change in any major way our predictions on the role 
of electrostatics. 

The bilayer consists of two undeformable and parallel 
leaflets, lying in the xy-plane and in contact, respectively, 
with the two monovalent salt solutions, as is shown in 
Fig. [TJ Let us label all quantities residing on the inner 
leaflet by the subscript 1 and those on the outer one by 
2. The inner leaflet (leaflet 1) is located at z = d and 
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FIG. 1. Schematic representation of the electrostatic potential 
for a binary lipid bilayer consisting of negatively charged and neu- 
tral lipids. The head groups of leaflets 1 and 2 are located at the 
plane z — dand z — 0, respectively. The (dimensionless) electro- 
static potentials in solutions 1 [z > d) and 2 (2: < 0) are denoted 
by •4'i{z) and tp2{z), respectively. The dimensionless surface po- 
tentials at the two leaflets are denoted by ^'i — ■il)i{z — d) and 
$2 = ^2(2 = 0), while ip-i^[z) is the (dimensionless) electric po- 
tential inside the hydrocarbon core region (0 < 2; < d) of the 
lipid membrane. The dielectric constant of the aqueous region 
(water) and lipid hydrocarbon core are denoted by and e-^, 
respectively. 

is in contact with the inner reservoir located at z > d, 
while the outer leaflet (leaflet 2) is located at z = and 
is in contact with the outer reservoir at 2: < 0. In some 
experiments spontaneous budding of charged domains to- 
ward the vesicular interior was reported . This can be 
explained by considering the increase in the osmotic pres- 
sure coupled with the change in spontaneous curvature 
of the outer leaflet. However, in the present study we as- 
sume that the bilayer remains flat and neglect the effect 
of the osmotic pressure generated across the bilayer. Our 
model is further based on three principal assumptions: (i) 
Because we do not consider membrane undulations and 
their curvature, the two leaflets have the same area A. 
(ii) The neutral and negatively charged lipids have the 
same cross-sectional area per lipid S. This is supported 
by many experiments on fluid-like membranes, (iii) Fur- 
thermore, each leaflet is taken as an incompressible two- 
dimensional fluid. Hence, the two leaflets consist of the 
same total number of lipid molecules N = A/Tj. 
The total free energy per lipid molecule is given by 

ftot{(t>l,4>2) = /mix(0l) + /mix('/'2) + /cl(01, ^2), (1) 

where all energies are measured in units of k^T (fee is 
the Boltzmann constant and T the temperature) , and (pi 
and (/)2 are the mole fractions of the negatively charged 
lipid in leaflets 1 and 2, respectively. The free energy 
/tot consists of two non-electrostatic /mix terms (one for 
each leaflet), and an electrostatic one /d that induces an 
inter-leaflet coupling. 
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The non-electrostatic /mix is the Flory-Huggins free 
energy of lateral mixing in each of the leaflets separately 



fmixi(f>) = (/"Inf/) + (1 ■ 



iln(l 



+ X0(l-0), (2) 



where x is the (non-electrostatic) interaction parameter 
between the two different lipids. Quite generally, x is 
taken to vary inversely with the temperature. The free 
energy /mix is a sum of the ideal entropy and enthalpy of 
mixing between the two lipid species. In the absence of 
the electrostatic interaction, fd — 0, the two-phase coex- 
istence region in the (x,</')-plane is delimited by a demix- 
ing curve that terminates at the critical point: Xc — 2 
and 4>c = 0.5. 

The electrostatic free energy /ei can be calculated 
through the charging process [23| : 



/c: 



02) = 



«'i(si,0) dsi 



, S2)dS2 



(3) 

where = e$i(z=d)/fcBT and *2 = e$2(z=0)/fcBT 
are the dimensionless surface potentials on leaflets 1 
{z = d) and 2 {z = 0), respectively, and e is the ele- 
mentary charge, while the electric potential inside the 
lipid hydrocarbon core region (0 < z < d) is denoted by 
<i>L(z). Within mean-field theory, the electric potential 
<i>(z) satisfies the Poisson-Boltzmann equation 
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where n is the salt concentration in bulk and the di- 
electric constant of the aqueous solution (water). Us- 
ing the dimensionless electrostatic potential = 
e^{z)/kBT, we obtain the following Poisson-Boltzmann 
equations in regions 1 and 2 for 'ipi^2iz), and Laplace 
equation in the core- region for ipi^i^), 

— Ki smh?/>i 



dz2 
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dz2 
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'. sinh ?/'2 



dz2 







z < 0, 



< z < d. 
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where n-^ ^ and K2 ^ 



and 



are the Debye screening length in 
2, respectively, defined by (^1,2)"^ = 



regions 1 

Y^£„fcBT/2e2ni,2- Using the Gauss theorem, the bound 
ary conditions of the two leaflets are given by 
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where S is the cross-sectional area per lipid head group 
and is the dielectric constant of the lipid hydrocarbon 
core-region. Far from both sides of the membrane, the 
bulk electric potential is taken to be zero, i/'iC-^ ^ 00) — 
and ^^2(2 — >■ —00) = 0. 



Under these boundary conditions, the above equations 
can be solved analytically, yielding transcendental ex- 
pressions for the two surface potentials, and '^2 



5'i = -2sinh ^ {pi(f)i + piH 
-^2 = -2sinh"^ {p2(t)2~P2HM!) 



(7) 



where A^* = ^'i — ^'2 is the membrane potential gap. 
Throughout the paper we will make use of three dimen- 
sionless quantities defined as: 



Pi = 

P2 = 

H = 



27r^ 
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where — / (Ans^k-BT) ~ tA is the Bjerrum length. 
The ratio P2/P1 — 1^1/ K2 is related to the ratio between 
the two Debye screening lengths, (^1.2)^^, of the two 
solutions, and H is the rescaled electrostatic coupling 
parameter between the two leaflets. By further defining 
(j>i = (pi + A^f and = '/'2 ~ A^ as effective mole 
fractions of the negatively charged lipid in each leaflet, 
we can rewrite Eq. ([7]) as: 



*i = -2 sinh ^(picpi) 
^2 = -2sinh"\p2 02) 



(9) 



Furthermore, the potential difference across the mem- 
brane A^* = v&i — ^1*2 satisfies a transcendental equation 



A^' = -2sinh ^(pi^^) + 2 sinh ^(p202 



(10) 



since (pi 2 themselves depend on A^P as defined above. 

By substituting back the surface potentials, Eq. ([7]), 
into Eq. ([3]), we obtain the final expression of the elec- 
trostatic free energy 
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The first term corresponds to the charging energy stored 
in the two charged leafiets that are analogous to a two- 
plate capacitor. The second and third terms represent 
the electrostatic energy of the two isolated charged mono- 
layers. In the absence of the electrostatic coupling be- 
tween the two leafiets {H = 0), /d is simply the sum of 
the free energies of two isolated leaflets, i.e., fci{(f>i, (f>2) = 
/m(0i) + fm{4>2)- Within the Poisson-Boltzmann theory, 
/m was derived in Ref. [2l|, and we briefly repeat this 



derivation in the Appendix. The result is 
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p = Pi for (f) — (f>i, and p = P2 for (f> = (f)2- We first solve 
Eq. ([To]) to obtain A^*, and then use it to calculate the 
total free energy from Eqs. (P), ©, (HB) and (fT2|) . 

In order to analyze the membrane stability toward a 
lateral phase separation, we calculate the spinodal sur- 
face given in the (0i, 02, x) parameter space by the con- 
dition 



should be equal for the two phases, i.e, = ffB- These 
conditions are sufficient to calculate the two-phase coex- 
istence region as well as the corresponding tie-lines in- 
side the coexistence region, as will be discussed below in 
Sec. mil We remark that the same conditions were also 
used in previous works [13, II3| ■ 



aVtot 



d'fu 



0. 



d(t)l J \ J \d(f)id<f)2, 
Using the total free-energy of Eq. ([TJ , we obtain 
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The first and second derivatives of the electrostatic free- 
energy, /oi(0i, 02), are given by Eq. (|TT|) and can be ex- 
pressed as 
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= 2 sinh ^(pi0i) = -^-i, 
2 sinh" ^(p2 02) = -^2, 



(15) 



and 
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where /;:;(0) = dVm(0)/d0 = 2p/^(p0)2 + 1. Finally, 
the spinodal surface Xsp(0i, 02) for fixed pi,p2,H can be 
determined according to Eqs. (fTO)) . (fT4l) and (IT6l) . 

The chemical potential is a function of the lipid con- 
centration 
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+ X(l-20i), 

+ x(l- 202). 

(17) 



Under the assumptions mentioned in the beginning in 
this section, the condition for phase coexistence (the bin- 
odal line) requires the equality of the chemical potential 
of the two phases, A and B, on each of the two leafiets, 
meaning that both fii^A = /^i,b and /i2,A = M2,b should 
be simultaneously satisfied on the two leaflets. In addi- 
tion, the thermodynamic potential 
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2) — /t( 



62) - Ml'/'l - M202 (18) 



A. The symmetric in/out reservoir case, pi = p2 

Let us first discuss the case of equal ionic strength 
{pi — P2) of the two solutions and review some of the 
results reported already in Ref. [11]. The critical point, 
because of the in and out symmetry, is located on the 
01 = 02 line, and the expression for the spinodal line is 
given by 

= 20(1 - 0) + Vl+PV+4pi?' ^^^^ 

where p — pi — P2 and = 0i =02- The minimum 
value of Xsp corresponds to the critical point. 

Let us take now explicitly the limits of high and low 
salt concentrations for the symmetric p case. The high- 
salt limit {p ^ 1) corresponds to the Debye-Hiickel 
regime for which the spinodal line is given by 

1 P 
^ 20(1 - 0) + 1 + ipH ' ^^^^ 

and the critical point is located at Xc = 2, 0^ = 0.5. 
This means that the bilayer behaves as if it were a neu- 
tral membrane when the salt concentration is sufficiently 
high. 

In the other limit of low salt (p ^ 1), the spinodal line 
is given by 

^^P"*" 20(1-0) + 0TIff' ^^^^ 

and it is independent of p. If the electrostatic coupling 
is weak enough {H <C 1), the respective critical point is 
located at 

Xc-(2 + V3)(l-^), 



■i-V3 4i^(^/3-l) 



(22) 



"^ 2 3 

In the absence of the electrostatic coupling {H — 0), 
the critical point coincides with that in the previous 
works [11 ill [13: 

= 2 + \/3~ 3.73, 

= (3- V3)/2~ 0.634. (23) 



Xc 



III. GLOBAL PHASE DIAGRAMS AND 
CRITICAL BEHAVIOR 

Based on the model described in the previous section, 
we present results for the critical point and phase dia- 
grams in the case of two coupled and charged leafiets. In 
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FIG. 2. (a) The critical interaction Xc and (b) critical com- 
position <j)c as a function of p on a semi-log plot for different 
H values in the symmetric case (pi = p2 = p). The solid, 
dashed, and dot-dashed lines are for H = 0.005,0.01,0.05, re- 
spectively. For H = 0, the critical values are represented by 
the dotted lines located at (a) Xc = 2 -I- -\/3 ~ 3.73 and (b) 
(?!>c = (3-^)/2~ 0.634. 



particular, we focus on the difference between a symmet- 
ric case in which the two solutions have the same ionic 
strength (pi = P2) and the asymmetric case in which 
they differ (pi p-i)- 

We start by giving some estimates for the physical 
choices of the parameter values of pi, p2, and H. The 
cross-sectional area of a lipid headgroup is set to about 
E ~ 65 (taken to be the same for the two lipids) , and 
the membrane thickness is c? ~ 50 A. From these values 
together with e„ ~ 80eoj — 3eo (eo is the vacuum per- 
mittivity) and T ~ 300 K, we obtain p ~ 0.67k~^ where 
kT^ is measured in Angstroms, and H — 5.6 x 10~*. 
In the following, we use a range of values p — 0.5, 5, 50 
and extended the H range to cover H = 0.005, 0.01, 0.05. 
In the case of monovalent salt, the values p = 0.5, 5, 50 
correspond, respectively, to = 0.75, 7.5, 75 A or to 
n = 16 M, 160 mM, 1.6 mM. Since the electrostatic cou- 



pling parameter H is proportional to S-^/d^ the larger H 
values correspond to smaller d and/or larger e^. Obvi- 
ously, the range in the p parameter of 10^ gives a range 
of lO'* in ionic strength, which probably is too large to 
compare with experiments. Hence it should be taken just 
to indicate trends with changing the salinity. 

We show first the phase separation in asymmetric 
charged bilayers sandwiched between the two ionic so- 
lutions having the same ionic strength (pi — P2), repeat- 
ing the results of Ref. [ll| . Then, we present results for 
asymmetric charged leaflets that are in contact with two 
ionic solutions having different ionic strengths (pi 7^ ^2)- 

A. Symmetric in/out reservoirs: pi = P2 

In Fig. [21 we show the critical interaction and compo- 
sition (xcj 0c) as a function of p for the symmetric reser- 
voir case, p — pi = P2- Both Xc and (pc are obtained by 
numerically minimizing Eq. (|19p. The solid, dashed, dot- 
dashed lines correspond, respectively, to H = 0.005, 0.01, 
and 0.05. The horizontal dotted lines located at Xc = 
2 + ^ 3.73 and 0c = (3 - V3)/2 = 0.634 are the 
critical values for the iJ — > limit and for high-salt con- 
ditions. 

When x> Xc the membrane undergoes a lateral phase 
separation, whereas for x < Xc the membrane is in a sin- 
gle (homogeneous) phase. For all H values, we find that 
Xc becomes smaller as p decreases. Namely, the phase- 
separation temperature is increased and the two-phase 
region has a larger extent for smaller p. This is because 
the Coulombic repulsion between the charged lipids be- 
comes too weak to overpower the attractive interaction 
that drives the phase separation. When H is large, Xc be- 
comes much smaller, as will be discussed later in Sec. lIVI 
In the high salt limit or when p is small enough, the limit 
of Xc — > 2 and (pc ^ 0.5 can be seen, as was discussed 
above in Sec. [ill In the opposite low salt limit, Xc and 
0c approaches those of Eq. (j^ . 

B. Asymmetric in/out reservoirs: pi 7^ P2 

The calculated value of Xc for the asymmetric case, in 
which the two solutions have different ionic strengths, is 
plotted in Fig. [3] Since Xsp equals to Xc at the critical 
point, Xc is obtained by calculating numerically the mini- 
mum value of Xsp via Eq. . The horizontal axis is the 
ratio P2/P1 in logarithmic scale, where P2/P1 = 1 corre- 
sponds to equal ionic strength of the two reservoirs. The 
solid, dashed, and dot-dashed lines denote fixed pi — 0.5, 
5, and 50, respectively, while the coupling parameter H 
has three different values H — 0.005, 0.01, and 0.05 in the 
three figure parts. Similar to the findings of Fig. [21 Xc 
decreases as the coupling H increases. This dependence 
of Xc on H is similar as in the symmetric case and is fur- 
ther discussed in Sec. IIVI In all plotted cases there is a 
peak around P2/P1 = 1, without a break in its slope there 
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H = 0.005 H = 0.01 H = 0.05 




P2/P1 P2/P1 P2/P1 



FIG. 3. The critical value Xc as function of the asymmetry parameter, P2/P1, which also is equal to the ratio of the two Debye 
lengths, ki/k2- The solid, dashed, dot-dashed lines are for pi — 0.5,5,50, respectively. The values of the H parameter are: (a) 
H = 0.005, (b) H = 0.01, and (c) H = 0.05. 



even for the solid line. For P2/P1 < 1, Xc decreases as 
the salt concentration in solution 2 is increased, similar 
to the symmetric case. 

Moreover, the Coulombic repulsion between charged 
lipids in the same leaflet becomes weak and Xc decreases 
when P2/P1 becomes smaller. For P2/P1 > 1, Xc is found 
to decrease and its value is more pronounced when H 
becomes larger, although the electrostatic screening is 
weaker. The physical origin of this behavior will be fur- 
ther discussed in Sec. ITVl 

The calculated phase diagrams in the (0i,02) plane 
for X = 3.8 and 4.0 are shown in Figs. [Hand [3 respec- 
tively. We plot all phase diagrams for a fixed pi — 5, 
while varying p2 and H . As can be seen in Fig. [Sj the 
system undergoes a phase separation for x values larger 
than Xc- The solid lines are the tie- lines within the co- 
existence region, while the dashed lines correspond to 
the spinodal lines. Each pair of points connected by the 
tie- line satisfies the conditions /ii.A = Mi,b, M2.A = M2,b 
[Eq. (fT7| ] and = 5b [Eq- (flS]) ] while the spinodal lines 
are obtained from Eqs. ([TU]). ([TUl and p^ . In some of 
the figures, inner and external spinodal lines can be seen. 
However, the inner lines are preempted by the external 
ones. Since Fig. IDJb), (e) and Fig. [5l[b), (e) are for the 
symmetric case {pi = P2), these four phase diagrams are 
symmetric with respect to the diagonal line, (f)i = 02- 
All other phase diagrams correspond to asymmetric sit- 
uations, pi ^ P2- 

The phase-coexistence region becomes relatively larger 
when P2/P1 ^ 1, corresponding to strongly screened sys- 
tems. However, Xc for the asymmetric case is smaller 
than that of the symmetric case. When H — 0.005, most 
of the tie-lines are nearly parallel to either 0i- or 02-axis. 
This implies that the phase separation in each leaflet 
takes place almost independently and without any no- 
ticeable correlation with the second leaflet. On the other 



hand, for H = 0.05, the tie-lines are tilted, indicating a 
strong coupling between the two leaflets. 

The three-phase-coexistence regions are indicated by 
shaded triangles in the phase diagrams. The tie-lines 
in the vicinity of the three-phase-coexistence regions are 
almost orthogonal to the principal diagonal = (t>2, in- 
dicating that 01 is small when (f)2 is large and vice versa. 
This occurs when the local concentration of charged lipid 
in one of the leaflets is increased, while the concentration 
in the second leaflet is decreased conversely. This behav- 
ior is caused by the electrostatic coupling between the 
leaflets, as is discussed in the next section. 

IV. DISCUSSION AND COMPARISON WITH 
EXPERIMENTS 

Several points merit further discussion. First, we com- 
pare our theoretical results with experiments, where it 
was reported that the phase separation is enhanced when 
salt is added to charged bilayers [l^. Such a tendency 
is well reproduced in Figs. 2] and [S] when the salt con- 
centration in solution 2 is increased or equivalently p2 is 
decreased. A similar trend was also presented in previ- 
ous works [T3-[l^ for charged membranes without any 
electrostatic inter-leaflet coupling. In addition, it was 
reported by Vequi-Suplicy et al. (Tlj that the phase- 
separation temperature increases in the presence of salt. 

In previous experimental studies p^ . [ll[ , vesicles were 
first dispersed in the aqueous solution, and only then 
salt was added, affecting only the ionic strength of the 
solution outside the vesicles. For this experimental pro- 
cedure the value of P2/P1 is smaller than unity (recalling 
that p2 corresponds to the outer reservoir). This situ- 
ation is presented in Fig. [3] where Xc ^ l/^c decreases 
as P2IP1 < 1 becomes smaller, in accord with the ex- 
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p, = 0.5 



P2 =50 




FIG. 4. Phase diagrams in the (0i,</>2) plane for x = 3.8 and pi = 5. The solid lines are the tie-lines and the dashed lines are 
the spinodal lines. The three-phase-coexistence regions are denoted by shaded triangles. The values of the other parameters are: (a) 
H = 0.005 and p2 = 0.5; (b) H = 0.005 and p2 = 5; (c) H = 0.005 and pa = 50; (d) H = 0.05 and p2 = 0.5; (e) H = 0.05 and 
P2 = 5; and (f) H = 0.05 and pa = 50. 



perimental observations. In addition, we find that also 
for P2/P1 > 1, Xc is decreased even if the electrostatic 
screening is weak. In order to confirm this behavior in 
experiments, only the outer solution should be diluted in 
order to yield P2/P1 > 1. Controlled experiments done 
on asymmetric charged membranes with asymmetric in 
and out ionic strengths are needed to confirm the predic- 
tions of the present work. Because the lipid composition 
in each leaflet prepared by gentle hydration or electrofor- 
mation is not well controlled, it is better to prepare vesi- 
cles by transferring water-in-oil droplets coated by lipids 
from an oil phase to a water phase [2^ HI] or by using 
asymmetric Montal-Mueller planar bilayers [H, [2^ . 

In the following, we discuss the physical reason why Xc 
decreases for P2/P1 > 1- Such a tendency is more pro- 
nounced for larger H ~ El/^, because the electrostatic 
coupling between the two leaflets becomes stronger. No- 
tice that the electrostatic coupling depends on pi and p2 
as well as on H. As schematically presented in Fig. [HI the 
electrostatic coupling is weaker when the salt concentra- 
tion in solution 2 is higher (upper panel of Fig. and 
the phase separation is suppressed. 

On the other hand, when the salt concentration in so- 
lution 2 is low (lower panel of Fig. [5]), the screening is 



not as efficient, and phase separation is induced due to 
the strong electrostatic coupling. Note that in the latter 
case, the strong demixing in leaflet 1 is induced by chang- 
ing the salt concentration in other reservoir (solution 2), 
which is coupled directly with the opposite leaflet. This is 
indeed an interesting situation because the higher salinity 
in solution 1 relative to 2 triggers the phase separation in 
leaflet 1 due to the strong inter-leaflet electrostatic cou- 
pling. The region where the tie-lines are almost parallel 
to the (jji ^xis means that the phase separation mainly 
takes place in leaflet 1. As can be seen in Fig. 0] (c) and 
(f) as well as in Fig.[5](c) and (f), such an region becomes 
much larger than the region where the phase separation 
strongly occur in leaflet 2 for large p2- On the other 
hand, the phase separation also occurs in leaflet 2, al- 
though the concentration difference between two phases 
is very small. As shown in Fig. [31 Xc is slightly decreased 
when pi is smaller, and can be understood by noting 
that the screening sufficiently weakens the electrostatic 
coupling. 

In Figs. [2] and [3] we have shown that Xc is smaller for 
larger H, i.e., when the electrostatic interaction across 
the membrane is large. Because of the electrostatic cou- 
pling, the local charge accumulation in one of the two 
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FIG. 5. Phase diagrams in tiie (0i,(/)2) plane for x = 4.0 and pi = 5. The solid lines are the tie-lines and the dashed lines are 
the spinodal lines. The three-phase-coexistence regions are denoted by shaded triangles. The values of the other parameters are: (a) 
H = 0.005 and p2 = 0.5; (b) H = 0.005 and p2 = 5; (c) H = 0.005 and p2 = 50; (d) H = 0.05 and p2 = 0.5; (e) H = 0.05 and 
P2 = 5; and (f) H = 0.05 and p2 = 50. 



leaflets suppresses the charge accumulation in the other 
leaflet, and it is energetically unfavorable for the charged 
domains on each leaflet to face each other. As a result, 
a phase separation is induced in the other leaflet driven 
by the electrostatic coupling. 

Next we elaborate on the physical meaning of the slope 
of the tie-lines in Figs. |4] and El When the tie-lines are 
located close to the three-phase region (triangle) and/or 
when H is large, the tie-lines tend to be almost orthog- 
onal to the principal diagonal. These tie-lines connect 
two points on the binodal line that have a large compo- 
sitional asymmetry between 0i and 02. It indicates that 
the charged domains in one leaflet do not prefer to face 
the charged domains in the other leaflet. The model by 
Baciu and May [T^ as well as the present model deal 
with the compositional coupling between the two leaflets 
arising exclusively from electrostatic interactions. Con- 
sequently, it is energetically favorable for the charged do- 
mains on each leaflet not to face each other. 

In order to see the contribution of the electrostatic 
coupling clearly, we consider the electrostatic coupling in 
the symmetric case (pi = P2)- The free energy of the 
inter-leaflet interaction /coup is obtained by subtracting 
the contributions of individual leaflets from the total free 



energy 

/coup('/'l, </'2) = /tot(</>l, (/>2) - /mix('/'l) " /mix(<^2) 



fn 



(24) 



For small composition differences, the potential difference 
across the membrane is approximately given by 



A* 



2p 



4Hp + y/1 + p^cbl 



^2), (25) 



where (/)av — {4>i + 02)/2 is the average composition of 
the charged membrane. By substituting Eq. (|25p into 
Eq. (P^. the free energy of the inter- leaflet interaction 
can be expanded [1^ in terms of (f)i — 4>2, 



2Hp^ 



coup V 



1 + F 



,2^,2 , 
av 



■iHp^Jl + p'^(j)': 



(26) 

From Eq. (1^ . it is clear that /coup is negative and 
becomes smaller when the difference between and 02 
is large, because the coefficient of (0i — ^2)^ is negative. 
This is consistent with the diagonal tie-lines in the phase 
diagrams, and a similar tendency can be seen even for 
the more general asymmetric systems. 
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FIG. 6. Schematic illustration of the phase separation for an 
asymmetric charged lipid bilayer when p2/pi > 1. On the up- 
per panel, the salt concentration in solution 2 is high (but still 
obeying P2/P1 > !)■ On the lower panel, the salinity in solution 
2 is reduced resulting in an increase of p2/pi- The electrostatic 
coupling becomes then stronger and the phase separation is in- 
duced. 



In general, the physical origin of the coupling be- 
tween the two leaflets can vary and does not have to 
rely exclusively on electrostatics interaction. Examples 
of such non-electrostatic effects were considered in pre- 
vious works [2^, [13] and include cholesterol flip-flop, dy- 
namic chain interdigitation, van der Waals interaction or 
composition-curvature coupling. The coupling between 
the two leaflets can be expressed generally by 



/coup = A(0i - 02 )^ 



(27) 



as shown in Eq. (1261) . The coupling coefficient. A, be- 
comes negative in the case of the electrostatic coupling. 
However, Wagner et al. [l^ analyzed a model for coupled 
bilayers containing such a phenomenological term as in 
Eq. (P71) but with A > 0. It is, therefore, important to 
consider the sign and typical value of A for different sys- 
tems containing electrostatic and non-electrostatic cou- 
plings and to reveal how the interactions between the two 
leaflets contribute to the coupling constant, A. Further 
experimental and theoretical investigations are required 
to gain more insight on this coupling and its origin. 

The present work is concerned with macro-phase sep- 
aration in lipid membranes. Membranes consisting of 
neutral and charged lipids can form also modulated 
phases [l^ that exhibit characteristic periodic structures. 
Hirose et al. ^2^ considered the coupling between two 
micro-phase separated bilayer leaflets. Various patterns. 



phases and their phase transitions were predicted by as- 
suming a phenomenological coupling between the two 
leaflets. 

We would like to close by mentioning again the two 
main approximations used in the present model, (i) We 
have treated the two leaflets as undeformable flat sheets. 
By doing so we neglected any possible bilayer deforma- 
tion that can occur because of osmotic pressure differ- 
ence resulting from a difference in ionic concentrations 
(Pi 7^ P2) on the two membrane sides. In experiments, 
phase-separated domains produce budding toward the in- 
terior of the vesicle and this is a direct consequence of the 
osmotic pressure on the spontaneous curvature ^Oi]. In 
future studies it may be of advantage to consider the 
combined effect of the phase separation and membrane 
deformation, (ii) We have used the Poisson-Boltzmann 
theory for symmetric monovalent salt in order to describe 
the electrostatic potential. It is known that such a mean- 
fleld theory does not treat correctly ion-ion correlation, 
especially when multivalent cations are involved. This 
important effect should also be taken into account in 
follow-up modeling. 

In conclusion, we have considered the phase separation 
in a charged and asymmetric lipid bilayer, located be- 
tween two ionic solutions having different ionic strengths. 
Considering the electrostatic effects on the phase separa- 
tion, we have introduced the electrostatic coupling be- 
tween the two leaflets. We studied the effect of the ionic 
strength difference between the two solutions, pi ^ P2, 
and the electrostatic coupling, H, on the membrane 
phase behavior. When the electrostatic coupling H is 
introduced, Xc ^ l/^c becomes smaller than for the un- 
coupled case. Moreover, even if the screening effect is 
small, Xc decreases by increasing the asymmetry in the 
salt concentrations. An increase in the ionic strength 
leads to an enlarged coexistence region in the phase dia- 
gram. 
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Appendix: Surface potential and charging free 
energy of an isolated leaflet 

1. Leaflet electrostatic potential 

We present the derivation of the surface potentials ^1^2 
[Eqs. d?]) and ([9])] and the electrostatic free energy of 
an isolated leaflet [Eq. ([T2|) ]. In order to find VPi, we 
multiply the Poisson-Boltzmann equation (g]) by "01 (^) 
and integrate it once, yielding 



V'lV'i dz = Ki / t/j'i sinh?/;i dz 



(A.l) 



where ip[ — d'i/'i/dz and Vi — d^ipi/dz^. Using the 
appropriate boundary conditions results in the following 
expression: 



^p'l = —2ki sinh [ 



and ai z = d. 



-2ki sinh ( — — 



(A.2) 



(A.3) 



Using Eq. (O and i}[(d) = ■i/'2(0) = A^/d, it follows that 
snih — 



^^A*. (A.4) 

KiS KiS Ane^iB d 

From this expression, the surface potential ^I'l is 

*i = -2 sinh"^ {pi(f>i + piHA^) , (A.5) 

where pi — 2tt£b/ kiY., H = Yi£^/ Ane^lBd as well as 
other parameters were introduced in Sec. |TT1 The surface 
potential on the second leaflet, ^'2, can be obtained in 
the same way. 



From this expression, the surface potential ^ = i){z = 0) 
is given by 



* = -2sinh"^(p</«). 



(A., 



Through the charging process, the electrostatic free- 
energy of an isolated leaflet is obtained by the following 
integral over ^I^ yielding an analytical expression for /m: 



20 



'^{s)ds 



In [p(j) + v/WP+T 

(A.9) 



2. Electrostatic free-energy of an isolated leaflet 

Next we calculate the electrostatic free-energy of an 
isolated leaflet, Eq. (fT2)) . using the charging method. A 
similar derivation can be found, e.g., in Refs. [l^, [2ll |. 
This energy corresponds to the energy that a leaflet feels 
in solution without the electrostatic coupling to the sec- 
ond leaflet. When a charged leaflet at z = is in contact 
with an ionic solution situated at z > 0, the Poisson- 
Boltzmann equation gives 

V>' = -2Ksinh . (A.6) 

This equation is the same as Eq. (jA.2p . The nonlinear 
first-order differential equation can then be integrated 
analytically: 



0(z) = -2 In 



1 



exp(Kz)coth[i sinh (p0) 



(A.7) 
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